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Abstract 

We describe calculations of cc production to next-to-next-to-leading order (NNLO) 
and next-to-next-to-leading logarithm (NNLL) near threshold in pp and ir~p interactions. 
We study the relevance of these calculations for existing cc total cross section data by 
examining their sensitivity to partonic threshold kinematics, their convergence properties 
and scale dependence. 



1 Introduction 



Calculations of charm production are still not under solid theoretical control. A good un- 
derstanding of the charm cross section is of particular interest for heavy ion physics. Charm 
production is an important contribution to the dilepton continuum in heavy ion collisions. In 
S+U and Pb+Pb interactions at the CERN SPS, with = 19.4 and 17.3 GeV/nucleon pair 1 
respectively, the dilepton yield in the mass range 1.5 < m fM +^- < 2.5 GeV is enhanced by a 
factor of 2-3 over the extrapolated proton-nucleus, pA, yield [1 . These data have been inter- 
preted as an enhancement of the charm production cross section in the system created in the 
heavy ion collision. Another possible source of the dilepton enhancement is thermal dilepton 
production, predicted in a quark-gluon plasma j2]. Finally, the total charm rate would be a 
useful reference for J/tp production in heavy ion collisions |2lEllI]- 

Although many future heavy ion experiments will be at high collider energies, \fS > 130 
GeV, some of the current and future experiments, like those at the SPS, are in the near- 
threshold region. The NA60 experiment was specifically designed to distinguish between charm 
decays and Drell-Yan-like production of dileptons j3] to determine the cause of the apparent 
enhancement. It will take heavy ion data at v^S = 17.3 GeV and pA data at vS = 29.1 
GeV. A new facility is being built at the GSI |6 a that will measure charm near threshold with 
= 6.98 GeV. 

Because the charm quark mass is a few times Aqcd, it is generally treated as a heavy quark 
in perturbative QCD calculations. However, its relative lightness results in a rather strong 
dependence of the total cross section on mass and scale, with up to a factor of 100 between the 
lowest and highest next-to- leading order (NLO) results [7j. There is also a rather broad spread 
in the measured charm production cross section data at fixed target energies. Much of this 
uncertainty arises from low statistics in the early experiments, assumptions of how much of the 
total charm yield results in final-state D mesons, and how the measured data are extrapolated 
to full phase space. The more recent data have improved considerably with new detection 
techniques and higher statistics. 

Improvements in the calculation of the charm cross section are difficult at all energies, but 
are perhaps possible when the cc pair is produced close to threshold, as we now describe. Fac- 
torization properties of QCD separate cross sections into universal, nonperturbative parton 
densities and a perturbatively calculable hard scattering function, the partonic cross section. 
Remnants of long-distance dynamics in the hard scattering function can dominate corrections at 
higher orders near production threshold. These Sudakov corrections have the form of distribu- 
tions singular at partonic threshold. Threshold resummation techniques organize these singular 
distributions to all orders, presumably extending the reach of QCD into near-threshold pro- 
duction. The singular functions organized by resummation are plus distributions, of the form 
[In' x/x] + 1 where x denotes the 'distance' from partonic threshold. At next-to-leading log (NLL) 
accuracy and beyond, proper account must be taken of the color structure of the hard scattering 
[HI E] for each partonic subprocess. 

Resummed cross sections are useful as generating functions for approximate finite order 
corrections to the cross section when expanded in powers of the strong coupling constant a s , 
as we do in this paper. The resummed cross sections may also be evaluated numerically. The 
charm fixed-target data were first compared to a leading log (LL) resummed calculation of the 

1 Unless otherwise specified, all pA energies are per nucleon and nucleus-nucleus energies are per nucleon 
pair. 
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total cross section in Ref. [TO] . Because the ratio m/A 3 is quite small, the expansion parameter, 
a s , is not and the LL resummation began to fail at y/S ~ 20 GeV. A NLL resummed evaluation 
in Ref. [11] found significant threshold corrections, albeit with a reduction in scale dependence. 

In this paper we work at finite order, using our results of Refs. [121 ESI- We have calculated 
the double-differential heavy quark hadroproduction cross sections up to next-to-next-to-leading 
order (NNLO), 0(aj), and next-to-next-to-leading logarithm (NNLL), i.e. keeping powers of 
the singular functions as low as I = 2i — 1 at order 0(a\ +3 ) where i — 0, 1, . . . ^21 E]- We only 
discuss QQ production in the ij = qq and gg channels since qg scattering first appears at NLO. 

Our studies focus on the kinematics of the proposed GSI facility and the CERN SPS proton 
and ion fixed target programs. We first briefly describe our NNLO-NNLL calculations in the 
next section. In section 3 we discuss whether the ion beam energies will actually produce charm 
in the near-threshold region where our calculations are in principle applicable. We show results 
for several values of the charm quark mass, m = 1.2, 1.5 and 1.8 GeV and for scales fi = m 
and 2m. We compare our results for the NNLO-NNLL inclusive cc cross section to charm 
production data and to the NLO cross sections in the relevant energy regime, and judge their 
value. Finally, we summarize our results in section 4. 



2 Resummation 

In our approach, the distance from partonic threshold in the singular functions depends on 
how the cross section is calculated. We either integrate over the momentum of the unobserved 
heavy quark or antiquark and determine the one-particle inclusive (1PI) cross section for the 
detected quark or treat the Q and Q as a pair in the integration, in pair invariant mass (PIM) 
kinematics. In 1PI kinematics, 

p{P 1 ) + p(P 2 ) — Q{ Pl ) + X{p x ) , (1) 

where X denotes any hadronic final state containing the heavy antiquark and Q(pi) is the 
identified heavy quark. The reaction in Eq. (0) is dominated by the partonic reaction 

i{h)+j(k 2 ) — > Q(pi) + X[Q](p' 2 ) . (2) 

At LO or if AT[Q](p' 2 ) = Q(p 2 ), the reaction is at partonic threshold with Q momentum p 2 . At 
threshold the heavy quarks are not necessarily produced at rest but with equal and opposite 
momentum. The partonic Mandelstam invariants are 

s = (ki + k 2 ) 2 , ti = (k 2 - pi) 2 - m 2 , ui = (ki - pi) 2 - m 2 , s 4 = s + t x + u x (3) 

where the last, S4 = (p 2 ) 2 — m 2 , is the inelasticity of the partonic reaction. At threshold, S4 = 0. 
Thus the distance from threshold in 1PI kinematics is x = s^/m 2 and the cross sections are 
functions of t\ and U\. In PIM kinematics the pair is treated as a unit so that, on the partonic 
level, we have 

i(h)+j(k 2 ) — QQ(p')+X(k>). (4) 

The square of the heavy quark pair mass is p' 2 = M 2 . At partonic threshold, X(k') = 0, the 
three Mandelstam invariants are 

M 2 M 2 
s = M 2 , h = — — (1-Pm cos9), ui = — — (1 + Pm cosO) (5) 
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where Pm = \A — 4m 2 /M 2 and 9 is the scattering angle in the parton center of mass frame. 
Now the distance from threshold is x = 1 — M 2 /s = 1 — z where z — 1 at threshold. In PIM 
kinematics the cross sections are functions of M 2 and cos 9. 

The resummation is done in moment space by making a Laplace transformation with respect 
to x, the distance from threshold. Then the singular functions become linear combinations of 
ln fc N with k < I + 1 and N = Ne lE where 7^ is the Euler constant. The 1PI resummed double 
differential partonic cross section in moment space is 



s 2 ,; J , l J = Tri HijP exp 
ati aui [ 



m/N 
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/i 



X 



S^P exp 



m/JV 



^rg(a.(AiO) 



(6) 



exp (jEi(N u , (jt, fi R )J exp (Ej(N t , /i R )J exp <j 2 y — (7; (a B (//')) + Tj («« (mO) 



MA 



where iV u = N(—Ui/m 2 ), N t = N(—ti/m 2 ), and (P) P refer to (anti-)path ordering. To find 
the PIM result, transform t\ and U\ to M 2 and cos# using Eq. (J3J and let = N t = N. The 
cross section depends on the 'hard', if^, and 'soft', SV,-, functions which are Hermitian matrices 
in the space of color exchanges. The 'hard' part contains no singular functions. The 'soft' 
component contains the singular functions associated with non-collinear soft-gluon emission. 
The soft anomalous dimension matrix, T l g, is two dimensional for qq and three for gg. The 
universal Sudakov factors are in the exponents Ei, expanded as 



exp(Ei(N u ,[x,m)) 
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These exponents also contain the effects of the singular functions due to soft-collinear radiation. 
The coefficients C l ^ n \ as well as the detailed derivation of the resummed and finite-order cross 



sections, can be found in Ref. ^3]. The momentum space cross sections to NNLO-NNLL are 
obtained by gathering terms at 0(a^) and 0(af), inverting the Laplace transformation and 
matching the iV-independent terms in and to the exact O(af) results in Refs. [Tl^ Ho]. 

Any difference in the integrated cross sections due to kinematics choice arises from the 
ambiguity of the estimates. At leading order the threshold condition is exact and there is no 
difference between the total cross sections in the two kinematic schemes. However, beyond LO 
additional soft partons are produced and there is a difference. To simplify the argument, the 
total partonic cross section may be expressed in terms of dimensionless scaling functions 
that depend only on rj = s/4m 2 — 1 [T3] . 



aij^s^m 2 ,^ 2 ) 



„ II ,,\ 00 k / ,,2 

aM J2^a s ^)) k j: /<*>(„) In''" 



fc=0 



1=0 



(8) 



We have constructed LL, NLL, and NNLL approximations to f-j' 1 ^ in the qq and gg channels 
for k < 2, I < k. Exact results are known for k = 1 and can be derived using renormalization 
group methods for k = 2, I = 1,2 Our calculations use the exact LO and NLO cross 

sections with the approximate NNLO-NNLL corrections. 

The inclusive hadronic cross section is obtained by convoluting the inclusive partonic cross 
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sections with the parton luminosity 

1 1 

^•(r,/! 2 ) = t J dx 1 Jdx 2 5(x 1 x 2 -T)(f) i/hl (x 1 ,iJ, 2 )(i) j/h2 (x2,iJ 2 ), (9) 



o o 



where (pi/h{x, fi 2 ) is the density of partons of flavor % in hadron h carrying a fraction x of the 
initial hadron momentum, at factorization scale \x. Then 



' dr 



/At 
— $y(r,/i 2 ) cr i:j (TS,m 2 ,fi 2 ) (10) 

h3=q,q,9 im 2 /s T 
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/ dlog 10 r/ — — ln ( 10 ) ®ij(v,V 2 ) ^ij(v,m 2 ,fi 2 ) 

_• _ = _ J —oo 1 + 7? 



where 



^=A- i = t^- 1 - ( n ) 

4m 2 4m 2 

Our investigations in Ref. ^3] showed that the approximation should hold if the convolution 
of the parton densities is not very sensitive to the high rj region. 



3 Numerical results 

In this section, we first test the applicability of our near-threshold treatment by calculating the 
parton luminosities. We then compare our approximate NNLO-NNLL cross sections with the 
exact NLO cross sections and the pp and n~p total charm cross sections. We also discuss the 
convergence properties and scale dependence of our results. 

3.1 Parton luminosities 

We first check if the parton luminosity peaks at low enough values of 77 for our near-threshold 
calculations to be applicable. Figures [T]|21 show the parton luminosities for v^S* = 6.98, 17.3 
and 29.1 GeV respectively, corresponding to the top new GSI energy, the CERN SPS lead 
beam, and the SPS proton beam. The results are shown for three values of the charm quark 
mass at each energy: m = 1.2, 1.5 and 1.8 GeV. Using the lowest value of m gives relatively 
good agreement with pp data at NLO when \i = 2m O [T71 [TH]. This mass also represents a 
lower bound on the possible range of charm quark masses. The charm quark pole mass is most 
commonly chosen as 1.5 GeV. Finally, 1.8 GeV represents an upper bound on the charm mass. 

The GRV 98 HO parton densities jTHj are used in the luminosity calculations, shown with 
scales /1 = m and 2m. The scale is not decreased below m because the minimum scale in the 
parton densities is larger than m/2. Such calculations would thus be of little value. 

Our results are primarily shown for the GRV 98 HO parton densities since they are consistent 
with the most recent pion parton densities by Gliick, Reya and Schienbein j2D] , denoted here as 
GRS2. We will also compare our results with those using CTEQ5M [21] since these densities 
were used in our previous calculations for heavier quarks |13j . 
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Figure 1: The parton luminosity for pp interactions at v S = 6.98 GeV as a function of rj using 
the GRV 98 HO densities. The left-hand side gives the qq luminosity, the right-hand side the 
gg luminosity. From top to bottom, the charm quark mass is 1.2 GeV in (a) and (b), 1.5 GeV 
in (c) and (d), and 1.8 GeV in (e) and (f). The solid curves are with \i = m and the dashed, 
/1 = 2m. 



At v^S* = 6.98 GeV, Fig. ^ charm production is well within the threshold region. There 
is no luminosity at rj > 2 for m = 1.5 and 1.8 GeV while rj = 3 is the highest rj with nonzero 
luminosity for m = 1.2 GeV. It is interesting to note that at even this rather low energy, 
the gg luminosity is still the largest, nearly a factor of ten greater than the qq for \x = m. 
The low qq luminosity can be attributed to the steeply falling antiquark density at large x, 
n/y/S ~ 0.2 - 0.6, over the charm mass range considered here. Each step in m drops the 
luminosity by an order of magnitude while moving the peak to lower rj, from rj ~ 0.35 when 
m = 1.2 GeV to ~ 0.15 when m = 1.8 GeV. Increasing the scale decreases the gg luminosity 
by a factor of w 3. 

The luminosity is considerably higher at = 17.3 GeV, the Pb+Pb center of mass energy 
at the CERN SPS, as shown in Fig. |21 Now the peak in the luminosity is at higher rj but still 
at rj < 1, even for m — 1.2 GeV. Thus, charm production may still be considered as within 
the threshold region. The location of the peak in r] is similar to that for it production at the 
Tevatron [T^]. Therefore, our calculations should be applicable. Changing the scale does not 
have a strong effect on the luminosity at this higher energy. The relative gg to qq luminosity 
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Figure 2: The parton luminosity for pp interactions at v S = 17.3 GeV as a function of i] using 
the GRV 98 HO densities. The left-hand side gives the qq luminosity, the right-hand side the 
gg luminosity. From top to bottom, the charm quark mass is 1.2 GeV in (a) and (b), 1.5 GeV 
in (c) and (d), and 1.8 GeV in (e) and (f). The solid curves are with \i = m and the dashed, 
fi = 2m. 

is smaller than at the lower energy. The higher SPS proton center-of-mass energy, ^ = 29.1 
GeV, leads to a luminosity peak at somewhat larger r], but remains less than unity, see Fig. El 
Since some weight is given to the region rj ~ 10, especially for the lowest m considered, this 
energy is the upper limit at which our calculation is applied. 

While all the calculations shown in Figs. ^01 have been made for the GRV 98 HO densities, 
the results for the CTEQ5M densities are quite similar. The luminosities in n~p interactions 
do differ however. At = 6.98 GeV, the qq luminosity is greater than the gg luminosity by 
a factor of ~ 1.5 for m = 1.2 GeV, increasing to nearly a factor of 5 for m = 1.8 GeV. By 
= 17.3 GeV, the situation has changed and the gg channel again dominates but only by 
a factor of 1.1-1.5. The smallest difference corresponds to the largest charm mass. Finally, at 
VS = 29.1 GeV, the 7r p and pp luminosities are rather similar. 

The scaling functions that contribute to the partonic cross section have been studied ex- 
tensively in Ref. [T3j. Since the scaling functions are essentially independent of m, we do not 
show them again here. We turn instead to a comparison of our calculations with the pp and 
7i~p total charm cross section data. 
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Figure 3: The parton luminosity for pp interactions at a/S = 29.1 GeV clS £L function of rj using 
the GRV 98 HO densities. The left-hand side gives the qq luminosity, the right-hand side the 
gg luminosity. From top to bottom, the charm quark mass is 1.2 GeV in (a) and (b), 1.5 GeV 
in (c) and (d), and 1.8 GeV in (e) and (f). The solid curves are with \i = m and the dashed, 
/x = 2m. 



3.2 Comparison with total charm data 

Comparisons of the NLO cross sections to the available cc data have been made to obtain 
the 'best' agreement by eye with the data by varying the mass, m, and scale, /x, for several 
combinations of m, /1, and parton density (THl [T71 [TH]. The best agreement with data at 
NLO is with /x = m = 1.3 GeV for the GRV 98 densities and with /x = m = 1.4 GeV and 
m = fi/2 = 1.2 GeV for CTEQ5M ^H]- Thus the hadroproduction data tend to favor rather 
light charm masses. 

We now turn to our NNLO-NNLL results. Since the NLO cross section is known exactly, 
we add the 0(af) NNLL approximate contribution to the exact NLO cross section in the qq 
and gg channels. We also apply a damping factor, 1 / -y/1 + rj, to temper the influence of the 
approximate scaling functions at large r] where we are further from threshold and have less 
control over our approximations. 

We find the greatest difference in the kinematics schemes for /x = m. The effect of the 
NNLO-NNLL terms is reduced for higher scales because a s (m) > a s (2m). The running of a s is 
significant for charm because the mass is rather low. Thus the strong running of the coupling 
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constant for these values of \i ensures a large scale dependence for charm production. There 
is also a significant parton density dependence in the results because /i/A 3 is not large. Since 
A 3 = 0.38 GeV for CTEQ5M and 0.2475 GeV for GRV 98 HO, the NNLO-NNLL corrections 
will be larger for CTEQ5M than for GRV 98 HO. 

The most important contribution to the NNLO-NNLL charm cross section is fg 2 '°\ shown 
in Fig. 5 of Ref. ^3] for both 1PI and PIM kinematics. The differences between fgg' ^ in the two 
kinematics at large rj are considerable. The functional dependence on rj begins to diverge for 
77 > 0.1. In 1PI kinematics, fgg' ^ is relatively small and positive until r] > 2 when it begins to 
grow. On the other hand, in PIM kinematics fg 2 ^ becomes large and negative with increasing 
rj. This large negative contribution can sometimes lead to a negative total cross section for 
H — m when the 0(aj) contribution is larger than the NLO cross section. The decrease in 
a s with fi = 2m keeps the total PIM cross section positive even though the NNLO-NNLL 
correction remains negative. 

We note here that there is some arbitrariness in the functional form of the gg scaling 
functions since the expression 1 — 2tiUi/s 2 , used in the gg Born cross section in Ref. [13], is 
equivalent at threshold to (t\ + u1)/s 2 . Either may be used in the scaling functions within the 
accuracy of our approximations. However, different choices can lead to non-trivial numerical 
differences in the cross section. We use the former choice here but point out that the latter 
choice increases the cross sections in both kinematics and leads to less negative (or even positive) 
values for the PIM results. We also note that including a class of subleading logarithms beyond 
NNLL, as derived in Ref. ^2], considerably lessens the difference between the NNLO qq scaling 
functions in the two kinematics over a large rj region. For the gg channel, especially given the 
ambiguity in the functional form of the gg scaling functions described above, it is harder to draw 
firm conclusions. A full evaluation of all the subleading terms requires two-loop calculations 
[2*2"] and has not yet been done. However, the class of subleading logarithms arising from the 
inversion from moment to momentum space has been calculated exactly [T2*] . 

As is shown below, the NNLO-NNLL cross sections in the 1PI and PIM kinematics are very 
different. Since including the subleading logs seems to make the scaling functions more similar, 
especially in the qq channel, we are motivated to present results using the average of the two 
kinematics as well as the individual 1PI and PIM results in the remainder of the paper. 

Our results with the GRV 98 HO densities are shown in Fig. [3] for pp interactions with 
both scales and the three charm quark masses. The approximate NNLO-NNLL cross sections 
are compared to the exact NLO results. We plot the 1PI and PIM cross sections as well as 
their average. The available data at y/S < 30 GeV are included. In these figures, we are not 
attempting to fit the data but to show the effect of the 0(a*) correction. The NNLO 1PI and 
PIM results with \i = m diverge most strongly from the NLO calculation. Increasing the charm 
mass reduces these differences, particularly for // = 2m. The PIM correction at NNLO-NNLL 
is larger and negative, giving a negative total cross section for fi = m = 1.2 GeV. 

The same trends are shown in Fig. [3] for the GRV 98 HO-based calculations of the total 
cross section in 7i~p interactions. We use the recently updated pion parton densities, GRS2 
|20j . This evaluation has a somewhat lower gluon density than the previous GRV-7T set [23 . 

We now show our pp results with CTEQ5M in Fig. [U] Since A 3 is larger for the CTEQ 
densities, all the cross sections are somewhat greater than those calculated with GRV 98 HO. 
The larger A 3 results in a larger a s and hence larger higher-order corrections. The NNLO 1PI 
and PIM results with \x = m diverge most strongly from the NLO calculation. The PIM result 
for m = 1.2 GeV is already negative at v^S* ~ 12 GeV. 
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Figure 4: The total cc cross section in pp interactions as a function of v^S" using the GRV 98 HO 
densities. The left-hand side employs the scale \x = m, the right-hand side, /i = 2m. From top 
to bottom, the charm quark mass is 1.2 GeV in (a) and (b), 1.5 GeV in (c) and (d), and 1.8 GeV 
in (e) and (f). The solid curves are the exact NLO result, the dashed curves, the approximate 
1PI NNLO-NNLL result, the dot-dashed curves, the approximate PIM NNLO-NNLL result, 
and the dotted curves, the average of the 1PI and PIM NNLO-NNLL results. 



The calculations agree only moderately well with the data. The earliest data are all rather 
low statistics and are mostly measurements of single D mesons. How the cc pairs hadronize 
is a particularly important question for energies near threshold where some channels may be 
energetically disfavored. We follow Ref. [7] and assume that since a(D s )/a(D° + D + ) ~ 0.2 
and a(A c )/a(D° + D + ) ~ 0.3, it is possible to obtain the total cc cross section from cr(DD) 
by multiplying it by ~ 1.5. This assumption could have a strong energy dependence so that as 
many charm hadrons as possible should be measured at each energy to study hadronization. 
The fragmentation into DD, DA C , A C D, A C A C , etc. could be studied at the GSI facility if the 
experiments are able to reconstruct charm mesons and baryons. 

Some of the total cross section data are based on lepton measurements. Semileptonic decays 
do not allow the momentum of the primary D meson to be entirely reconstructed, adding an 
additional layer of experimental uncertainty. The CERN SPS results are based on muon spec- 
trometers which cannot unambiguously determine the identity of the primary hadron. Finally, 
some of the data are taken on nuclear targets and then extrapolated to pp assuming that the 
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Figure 5: The total cc cross section in ir~p interactions as a function of \fS using the GRV 98 
HO proton densities and the GRS2 pion densities. The left-hand side employs the scale \i = m, 
the right-hand side, \x = 2m. From top to bottom, the charm quark mass is 1.2 GeV in (a) and 
(b), 1.5 GeV in (c) and (d), and 1.8 GeV in (e) and (f). The solid curves are the exact NLO 
result, the dashed curves, the approximate 1PI NNLO-NNLL result, the dot-dashed curves, the 
approximate PIM NNLO-NNLL result, and the dotted curves, the average of the 1PI and PIM 
NNLO-NNLL results. 



cross sections scale linearly with A, supported by fixed target measurements of the A depen- 
dence of charm production [23] . Thus the placing of data on the plots is primarily to guide the 
eye. 

The total cross sections in pp interactions for the three energies we have discussed are given 
in Table ^ The scale dependence is not necessarily reduced at NNLO-NNLL relative to the 
exact NLO. Due to the complete dominance of the gg channel in pp interactions, the dependence 
on kinematics choice is large. Because fjjg' ^ grows with r\ at large r\ in 1PI kinematics, the 1PI 
cross section is always larger than the NLO exact result. On the other hand, the PIM result 
always underestimates the exact NLO calculation since fj?'® is large and negative at large rj 
in PIM kinematics. Thus the PIM cross section with [i = 2m can be greater than that with 
fi = m. The average of the 1PI and PIM cross sections typically remains somewhat above the 
NLO result. 

Recall that all the results include the damping factor. Without the damping factor, the 
0(af) contribution can increase up to 40% at the highest energies. 
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Figure 6: The total cc cross section in pp interactions as a function of \fS using the CTEQ5M 
densities. The left-hand side employs the scale \x = m, the right-hand side, /i = 2m. From top 
to bottom, the charm quark mass is 1.2 GeV in (a) and (b), 1.5 GeV in (c) and (d), and 1.8 GeV 
in (e) and (f). The solid curves are the exact NLO result, the dashed curves, the approximate 
1PI NNLO-NNLL result, the dot-dashed curves, the approximate PIM NNLO-NNLL result, 
and the dotted curves, the average of the 1PI and PIM NNLO-NNLL results. 

Note that even for the NLO exact cross section, the CTEQ5M results are larger than the 
GRV 98 HO results. Because A 3 is greater for CTEQ5M, the NNLO-NNLL corrections will be 
most important for this set. The dependence on parton densities is reduced for heavier quarks 
since the ratio fi/A ^> 1 and a change in the ratio will not affect the value of a s significantly for 
large /i. That is a primary reason why the scale dependence is small for top quark production 
|12[ I13|. In this case, the running of a s is primarily responsible for the dependence of the results 
on parton density and the large kinematics dependence. 

As discussed before, attempts to find agreement with the pp and ir~p charm production 
data by varying the mass and scale usually result in a rather low value of m at NLO, ~ 1.2 
GeV for fi = 2m. The NNLO-NNLL results, particularly with 1PI kinematics, suggest that the 
same general energy dependence may be obtained for higher masses and lower scales so that a 
value of m closer to the pole mass, 1.5 GeV, would agree with the fixed target data. Using the 
GRV 98 HO densities, the NLO calculation with m = 1.2 GeV and /i = 2m is in reasonable 
agreement with the approximate NNLO-NNLL 1PI result obtained with n = m = 1.5 GeV. 
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Even at NNLO-NNLL, m = 1.8 GeV still underestimates the data considerably. 

Finally, we briefly compare our finite order result to the NLL-resummed result for pp — > cc 
shown in Fig. 15 of Ref. [TT]. Note that the resummed cross section depends on the prescription 
used to obtain the momentum space result. Our NLO calculations with m = 1.5 GeV and 
CTEQ5M are in good agreement with their NLO calculation using MRSR2 [25} since these two 
sets have similar values of A. The NLL resummed cross section, calculated using the partonic 
total cross sections, is in rather good agreement with our kinematics-averaged NNLO-NNLL 
results. 

3.3 Convergence properties 

Given the large corrections at NNLO-NNLL, one can question if the calculated cross section will 
stabilize such that all data can be described with one value of the charm mass even if NNNLO 
and higher order corrections were calculated. Since the NNLO-NNLL 1PI results allow the data 
to be described relatively well with \i — m — 1.5 GeV instead of m = 1.2 GeV and fi = 2m, it 
is quite possible that further corrections could lead to agreement with the data for still higher 
values of m. We can at least partially address this issue through an investigation of the K 
factors. 

We first briefly discuss the calculation of the 'first order K factor', K^ = ctnlo/^lo- Since 
olo can be calculated with LO or NLO parton densities, the value of the K factor depends 
on which way <tlo is calculated. In Ref. K^ was defined in three ways. We compare 
two of the definitions here: Kq~\ where o"lo is calculated with NLO parton densities and a 
two-loop evaluation of a s , and K^\ where ctlo is calculated with LO parton densities and 
a one-loop evaluation of a s . In both cases, the Born and 0(a^) contributions to onlo are 
calculated with NLO parton densities and two-loop evaluations of a s . The first definition, 
Kq , indicates the convergence of terms in a fixed-order calculation while the second, K%\ 
indicates the convergence of the hadronic calculation towards a result. If Kq > K^\ the 
convergence of the hadronic cross section is more likely. 

In Fig. [7[ we compare K^p and for charm production in pp interactions with the GRV 
98 HO densities. The dotted curves are Kq and the solid curves are K£\ Here, > 
for all values of fi and m shown. We see that for /i = m, Krp does not change much with m. 
It is « 3.5 for y/S = 6.98 GeV, decreasing to « 1.9 at y/S = 29.1 GeV for all cases. There is 
slightly more change with m for /j, = 2m since K% tends to flatten at lower y/S with increasing 
mass. At higher y/S, K% is somewhat bigger for the larger scale, ~ 2.2. On the other hand, 
Kq is a stronger function of both \i and m. When \i — m, increases from w 3.7 for 
m = 1.2 GeV to 4.2 for m = 1.8 GeV at y/S = 6.98 GeV, decreasing to w 2.2 for higher y/S. 
When H = 2m, K { 1) increases 10-15%. 

If, instead, we use the CTEQ5M densities to calculate the NLO K factors in pp interactions, 
> at all energies and the values are all larger than for the GRV 98 HO densities. The 
larger K factors for CTEQ5M are due to the larger a s with the CTEQ5M A3. 

More interesting are the K factors in ir~p interactions, shown in Fig. |H| With the GRV 98 
HO densities, again > K2 for & H A* an d m but the values are much smaller and almost 
independent of energy. When // = m, K2 ~ 1.5 and Kq ~ 2.2 while when \x = 2m, K^ ~ 2 
and Kq ~ 2.5. Both the reduced size of the K factors as well as their near energy independence 
can be attributed to the greater importance of qq annihilation in charm production, particularly 
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Figure 7: The theoretical K factors for cc production in pp interactions as a function of y/S 
using the GRV 98 HO densities. The left-hand side employs the scale /i = m, the right-hand 
side, fi = 2m. From top to bottom, the charm quark mass is 1.2 GeV in (a) and (b), 1.5 GeV 
in (c) and (d), and 1.8 GeV in (e) and (f). The solid curves are K% \ the dotted curves, Kq , 
the dashed curves, K[p\ } the dot-dashed, Kpj U , and the triple dot-dashed, Kj$ e , as defined in 
the text. 



at low V~S. The O(a^) contribution to the qq channel is considerably smaller relative to ctlo 
than in the gg channel, reducing the K factor. Since the pp cross sections are dominated by 
the gg channel, their K factor is always large. 

All these NLO K factors are rather large for pp interactions, typically a factor of two or 
more. We can define the next order K factor, = onnlo-nnll/onlo, and see if it is reduced 
relative to K^. We define K[% where the NNLO-NNLL IPI cross section is in the numerator, 
Kpi M , with the PIM cross section in the numerator, and Kj^ c where the numerator is the 
average of the IPI and PIM cross sections. We cannot distinguish between calculations with 
parton densities at different orders for these K factors, as we did for K^ l \ because there are 
no full NNLO parton densities available. Now both numerator and denominator are calculated 
with all NLO parton densities and the two-loop evaluation of a s . By calculating for the 
NLO parton densities, we can identify improvements in the fixed-order results relative to Kq . 

Figure [7| also shows K^p\, K P i M and Kj® in the dashed, dot-dashed and triple dot-dashed 
curves respectively. The IPI K factor, k[ p \, is not strongly dependent on /i, m or VS. It 
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Figure 8: The theoretical K factors for cc production in n~p interactions as a function of 
\fS using the GRV 98 HO proton densities and the GRS2 pion densities. The left-hand side 
employs the scale fi = m, the right-hand side, \x = 2m. From top to bottom, the charm quark 
mass is 1.2 GeV in (a) and (b), 1.5 GeV in (c) and (d), and 1.8 GeV in (e) and (f). The solid 
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curves are ii 2 \ 



the dotted curves, Kq~\ the dashed curves, K\pj, the dot-dashed, Kpj M , and 



(2) 



'(2) 



the triple dot-dashed, K^) e , as defined in the text. 



decreases slowly from 3 at low v^S* to m 2 at higher \/~S. Consequently for /i — m, it is 
smaller than either calculation of at low v^S* and larger than as we move out of 
the threshold region. When \i = 2m, K^p\ is always less than However, Kpj M has the 

strongest /i, m and v^S dependence of all the K factors even though it is the smallest. The 
fact that it is smaller is less indicative of the ultimate convergence of the expansion than of 
the fact that the NNLO-NNLL PIM contribution to the cross section is large and negative. 
This is demonstrated by the change of Kp} M from positive to negative for > 17 GeV when 
H — m — 1.2 GeV. As the charm mass increases, the slope of Kpi M with decreases and 
the K factor remains positive. The average K factor, K^\, is less than both and 
everywhere. It remains positive and is not strongly dependent on m, yt, or \^S, remaining ~ 2—3 
at low \/S and decreasing to ~ 1 at higher v^S*- 

Similar trends are observed for the CTEQ5M densities but now the K^> factors are all 
larger and more dependent on ji and m, even Kfp\. Now also K$ M < for /j, — m — 1.2 GeV 
and 1.5 GeV at large y/S. For ir~p interactions, on the other hand, all the K factors are again 
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the K factors shown in Fig. [7[ we see that K\p( is ~ 10 at low \/~S and ~ 4 — 5 at larger \fS. 

K% at high y/~S. It is difficult to tell from these results if 



smaller and is almost independent of energy, as shown in Fig. |HJ 

We could form another K factor, K^ 2 '^ = onnlo-nnll/olo = K^K^ , to test the conver- 
gence of the hadronic cross section. The result is not complete because the NNLO-NNLL cross 
section is only approximate and the NNLO parton densities are unavailable. By multiplying 

Note that K^J is smaller and K^J ~ 
further, higher order, K factors such as onnnlo /cnnlo-nnll will be consistently smaller than 
2 near threshold or not. The convergence of the hadronic cross section at low scales is thus 
left in doubt although we note that subleading terms may have some effect on the convergence 
properties of the cross section [T2|. Since 

< for fi = 2m, one can expect that the 
next-order correction might be still smaller but this is not certain. In any case, it is not a 
guarantee that even larger masses may be needed to obtain the by-eye agreement at NLO with 
further higher-order corrections. 
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Figure 9: The scale dependence of cc production in pp interactions as a function of v^S*- The 
left-hand side employs the GRV 98 HO parton densities, the right-hand side, CTEQ5M. From 
top to bottom, the charm quark mass is 1.2 GeV in (a) and (b), 1.5 GeV in (c) and (d), and 1.8 
GeV in (e) and (f). The LO (dotted), NLO (solid), NNLO-NNLL 1PI (dashed), NNLO-NNLL 
PIM (dot-dashed) and the 1PI-PIM average (triple dot-dashed) ratios are shown. 
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Figure 10: The scale dependence of cc production in ir~p interactions as a function of v^S", 
using the GRV 98 HO proton densities and the GRS2 pion densities. From top to bottom, the 
charm quark mass is 1.2 GeV in (a) and (b), 1.5 GeV in (c) and (d), and 1.8 GeV in (e) and (f). 
The LO (dotted), NLO (solid), NNLO-NNLL 1PI (dashed), NNLO-NNLL PIM (dot-dashed) 
and the 1PI-PIM average (triple dot-dashed) ratios are shown. 



3.4 Scale dependence 

We now turn to a comparison of the scale dependence at LO, NLO and NNLO-NNLL, shown 
in Fig. El for the ratio a(fi = m)/a{^ = 2m) as a function of v^S in pp interactions. The GRV 
98 HO scale dependence is shown on the left-hand side and the CTEQ5M scale dependence 
on the right-hand side. The GRV 98 HO scale dependence is largest for the LO ratio. The 
NNLO-NNLL 1PI ratio is somewhat larger than the NLO ratio for m = 1.2 GeV but the two 
are almost identical for m — 1.8 GeV. Thus for 1PI, the scale dependence is not reduced relative 
to NLO but it is not really increased either. The NNLO-NNLL PIM ratio is smaller than the 
other ratios but it is a stronger function of energy, again becoming negative for m = 1.2 GeV 
as \/S increases. This ratio has the largest mass dependence. The ratio for the average of the 
NNLO-NNLL 1PI and PIM cross sections remains below the NLO ratio for all energies shown. 

The CTEQ5M scale dependence is considerably larger. In this case, when m = 1.2 GeV, 
the NNLO-NNLL 1PI scale dependence is larger even than the LO, dropping below it only 
when m = 1.8 GeV. The overall increase in the scale dependence relative to GRV 98 HO is 
again related to the larger a s for low scales with the CTEQ densities which also increases the 
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NNLO-NNLL contribution to the total cross section. 

The scale dependence of n~p production, shown in Fig. [TUJ is considerably reduced and, 
except for the NNLO-NNLL PIM ratio, nearly independent of energy. This reduction in the 
scale dependence can again be attributed to the relatively larger qq contribution to the total 
cross section. 



4 Conclusions 

We have studied the behavior of NNLO-NNLL calculations for charm production near threshold 
in pp and 7i~p interactions in both 1PI and PIM kinematics. We find that there are large 
differences in the results for different kinematics. Thus, the uncertainties in the cross sections 
remain large even at NNLO-NNLL. There are additional uncertainties, mentioned earlier, due 
to the functional form of the scaling functions and to subleading logarithms which have some 
effect on the values and convergence of the 1PI and PIM cross sections. We note that the 
average of the 1PI and PIM results has some nice properties such as better convergence and 
less scale dependence. For both pp and n~p interactions, either choice of scale, and all values of 
y/S (except when fj, = m = 1.2 GeV at high y/S) , the average of the NNLO-NNLL 1PI and PIM 
results lies above the NLO cross section. Thus, the charm mass need not be too low to agree 
with the data. However the poor convergence properties as well as other uncertainties of the 1PI 
and PIM results make any quantitative statement about the inclusive charm hadroproduction 
cross section difficult. 
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cx(NLO) (fib) tr(lPI) (/lb) 


(t(PIM) (/lb) 


(GeV) 


jj, = m 


jj, = 2m fj, — m 


/j, = 2m 


(i = m 


li = 2m 


m = 1.2 GeV 


GRV 98 HO 


6.98 


0.45 


0.11 1.29 


0.29 


0.70 


0.21 


17.3 


11.1 


5.61 26.5 


12.4 


0.075 


4.71 


29.1 


28.8 


18.8 69.0 


40.3 


-12.9 


8.09 


CTEQ5M 


6.98 


0.76 


0.13 3.50 


0.44 


1.48 


0.29 


17.3 


19.5 


7.33 75.5 


20.5 


-17.0 


6.32 


29.1 


51.9 


26.2 200.1 


70.7 


-101.4 


6.66 


m = 1.5 GeV 


GRV 98 HO 


6.98 


0.028 


0.0062 0.079 


0.017 


0.056 


0.014 


17.3 


2.75 


1.31 6.13 


2.83 


1.44 


1.44 


29.1 


9.44 


5.77 20.4 


11.9 


0.47 


4.00 


CTEQ5M 


6.98 


0.043 


0.0077 0.17 


0.025 


0.11 


0.019 


17.3 


4.13 


1.55 13.1 


4.08 


1.19 


1.86 


29.1 


15.1 


7.37 46.3 


18.5 


-8.50 


4.82 


m = 1.8 GeV 


GRV 98 HO 


6.98 


0.0014 


0.00030 0.0042 


0.00086 


0.0034 


0.00074 


17.3 


0.76 


0.35 1.63 


0.74 


0.63 


0.44 


29.1 


3.45 


2.13 7.05 


4.08 


1.24 


1.77 


CTEQ5M 


6.98 


0.0024 


0.00042 0.0094 


0.0014 


0.0070 


0.0012 


17.3 


1.01 


0.38 2.89 


0.97 


0.82 


0.54 


29.1 


4.97 


2.42 13.5 


5.79 


0.32 


2.18 



Table 1: Charm total cross sections in pp interactions at the three fixed target energies consid- 
ered. The results are given for the NLO exact, the 1PI and PIM NNLO-NNLL approximate 
cross sections. The cc cross sections using the GRV 98 HO and CTEQ5M parton densities are 
compared. 
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